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The Influence of Percolation in the generalized Chalker— Coddington 

Model 
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(Received ) 

We numerically investigate the influence of classical percolation on the quantum Hall 
localization-delocalization transition. This is accomplished within the framework of the gen- 
eralized Chalker-Coddington network model which allows us to control the number of classical 
saddle points by setting the width W of the saddle point distribution. It is found that increasing 
this width causes a new microscopic length scale to appear which depends on W and scales with 
the exponent X w 1.36 which indicates a close connection to the classical percolation length 
£ and its exponent v v = 4/3. Furthermore, the influence of an increase in W on the spectral 
statistics of the quasienergies of the network model is investigated. An effect similar to the 
increase of the potential correlation length in the Landau model is seen. 

KEYWORDS: Chalker-Coddington model, level spacing distribution, scaling, quantum Hall systems, multifractality 
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§1. Introduction 

The Chalker-Coddington networkB is a model for quan- 
tum Hall systems with long ranged disorder potentials. 
It represents a system of two-dimensional (2D) electrons 
in a strong magnetic field and smooth disorder poten- 
tial. As a member of the quantum Hall universality 
classEP it has been used to determine various critical 
quantities at the localization-delocalization (LD) tran- 
sition point bfitmeen the quantized plateaus of the Hall 
conductance.EfHlaB 

The model is based on the semi-classical time evolution 
picture of 2D electrons moving along the equipotential 
contours of a smooth disorder potential under the influ- 
ence of a strong magnetic field. The electronic states are 
defined by the amplitudes on the network links repre- 
senting the equipotential contours and the time evolu- 
tion is determined by unitary scattering matrices at the 
nodes of the network corresponding to the saddle points 
of the random potential where tunneling between con- 
tours occurs. The scattering strength is determined by 
the electron energy and the energy of the saddle point. 
In contrast to the original network model introduced by 
Chalker and Coddington, who explicitly excluded any 
percolation effects by setting all saddle point energies to 
zero, the generalized version of the model allows percola- 
tion effects by introducing an energy rajaeaJ— W, W] for 
the distribution of saddle point energies BtsB' While this 
generalization of the model does not change the critical 
behavior as long as the investigated systems are large 
enough, it introduces a microscopic length scale a that 
depends on W/Eu. where E t is the tunneling energy at 
the saddle pointsS* Since the motion of electrons with 
energy E at saddle points with energies Uk that obey 
\uk — E\ S> E t follows the classical path, it was concluded 
in refs. || and [?] that this length scale must be connected 
to the classical percolation length, meaning that a scales 



(1) 



with the classical percolation exponent v p = 4/3, i.e. 
Et) \E t 

As long as the system size L is much larger than this 
length scale the generalized model will show the same 
critical properties as the original. 

In the following we will first give a short description 
of the generalized Chalker-Coddington model followed 
by a review of the arguments given in refs. || and [?] for 
the influence of percolation effects. We will then give 
numerical evidence that the microscopic length scale in- 
deed shows the scaling behavior we expected. After that 
we take a look at the spectral properties of the network 
model and their dependence on W . We will show that 
the shape of the level spacing distribution function at 
criticality changes more and more to Poissonian behavior 
with increasing W. This is also the case for the number 
variance and consequently the spectral compressibility. 
In spite of these changes the scaling exponent v seems to 
be unaffected. The change of shape in the level spacing 
distribution is similar to that observed by Ono et al.cP 
when increasing the correlation length of the disorder 
potential. 

§2. The Network Model 

As we have already mentioned the Chalker-Coddington 
model is based on the semi-classical picture of electrons 
in two dimensions moving under the influence of a strong 
perpendicular magnetic field B in a long-ranged disorder 
potential V, i.e. the correlation length ly of V is large 
compared to the magnetic length l c — y^hc/eB. Start- 
ing from this picture the following network model was 
developed.ErlEP 

The network consists of a 2D regular lattice (see Fig. 1) 
whose links are unidirectional channels and whose nodes 
are scattering centers represented by 2 x 2 unitary scat- 
tering matrices S , where k is the node index. The ma- 
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Fig. 1. The Chalker-Coddington network model. At each saddle 
point the matrix S describes the scattering from incoming to 
outgoing channels. The network operator U performs the same 
task for the entire network. 

trices S k map the amplitudes of the two incoming links 
onto the two outgoing ones. The unidirectional links rep- 
resent the equipotential contours of the potential V and 
their random length is simulated by a random kinetic 
phase given to each link. The nodes represent the saddle 
points of V where the electrons can tunnel between con- 
tours. The tunneling amplitude T for each, saddle point 
is given by T — (1 + exp((E — Uk)/E t ))P where Uk is 
the energy of the saddle point with index k and E t is the 
tunneling energy. E t is of the form E t = if-c, where c 
is the average curvature of the random potential which 
can be approximated by c « V Q /ly with V = {V 2 ) 1 / 2 . 
Taking into account that the saddle point energies are 
randomly distributed in the interval [— W, W] gives .us 
the following approximation for the tunneling energy® 

E, = w(h-)\ (2) 

This relation connects the quotient E t /W which is the 
relevant parameter of the network model directly to the 
relevant parameter l c jly of other smooth disorder mod- 
els. It also connects the two limiting cases of the network 
model (W/Et — the original Chalker-Coddington case 
and W/E t = oo the classical network percolation case) to 
the uncorrelated potential (7y=0) and the classical mo- 
tion in a magnetic field and random potential (l c = 0) 
case. 

The system shows an LD phase transition when E ap- 
proaches the critical energy E c — (V) — 0. In this case 
the correlation length £ scales as £ cx \E — E c \ u , where 
v w 2.3. For large enough systems the correlation length 
exponent v does not depend on Wn 1 This is also true 
for other critical exponents like ao ~ 2.27 detennined by 
multifractal analysis of critical eigenfunctionsU 

A wave function or state of the network is defined by 
a normalized vector 'J = (tpi, ... ,ip n ) whose elements 
tpi are the complex amplitudes on the n = 2L 2 network 



links, where L is the system size. A wave function that at 
all nodes k of the network obeys the scattering condition 
ipl(i t k) = J2j Sijipi(j.k), where k) maps the matrix in- 
dex i £ {0, 1, 2, 3} of the scattering matrix S k at node k 
to the respective link of the network, is stationary under 
scattering and therefore an eigenfunction of the system. 

We can take all the matrix elements of the S k and 
arrange them into a single operator U , so that the sta- 
tionarity condition becomes 

U(E)V = (3) 

The network operator U also functions as a time evo- 
lution operator for the network states, i.e. "i>(t + x).-== 
U^(i), where r is the characteristic scattering timeoEP 
Equation (||) states that a wave function which is an 
eigenfunction of U(E) with eigenvalue 1 is an eigen- 
function of the modeled system. Such eigenfunctions 
will only occur at discrete values E n = E forming the 
eigenenergy spectrum of the system. These eigenener- 
gies are not easily determined, but it was found in rcf. [l2| 
that the eigenvalues uj a defined by the equation 

U(E)V a = e luJ «( E) y a (4) 

for a fixed value of E show the same statistics as the 
eigenenergies E n close to E. These so called quasiener- 
gies w a {E) are much easier to determine. This and the 
fact that one can choose the exact point on the energy 
scale where we want to investigate the statistics-make 
them the ideal tool for spectral investigations.E3E3 We 
can set E = E c and use all the critical quasienergies 
<jJa(E c ), i.e. the eigenphases of U(E C ) to determine spec- 
tral properties at criticality. This is an enormous ad- 
vantage compared to other methods where only a small 
fraction of the spectrum is critical. We can even improve 
this method if we use the fact that for every eigenphase 
uj a of the network operator the phase uj a + n (or in the 
case of double periodic boundary conditions the phases 
w a + 7r/2, uj a + 7r and oj a + 37r/2) is also an eigenphase 
(see Appendix ) . This leads to a reduction of the matrix 
size increasing the speed of the numerical determination 
and also the attainable system sizes. Additionally, ev- 
ery eigenstate determined for a unitary operator U (E c ) 
at the critical energy_E c is critical and can be used for 
multifractal analysisBeP 

§3. Percolation and Multifractality 

Although percolation, i.e. W > 0, does not seem to 
have an influence on scaling and thereby on the result 
of the multifractal analysis, there are nevertheless some 
effects which have to be investigated.© They are very 
obvious if one looks at the wave functions themselves 
(see Fig. 2), but also show some influence in the scaling 
analysis. 

The eigenstates in Fig. 2 obtained in the generalized 
network model at E — E c = show the typical self- 
similar shape of critical wave functions which lecLtp the 
use of multifractal analysis on critical systemsQ In 
this context the scale invariance of the eigenfunctions is 
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Fig. 2. Critical wave functions (E = 0) for W = and W = 150. 
The system size is 256x256 saddle points. Darker areas denote 

lower square amplitude. Fi S- 3 ' The scalin S of the avera § e box-probabilities for different 

values of W at q = —.5. System size L = 250. 



investigated by obtaining the following quantities: 

m{l h ) = [ dr|*(r + r )| 2 , (5) 

which denote the probability of a particle to be in a box 
of linear size fa centered at i"o , respectively. The disorder 
averaged g-moments of m(lb), m q (lb) — (m q (fa)), scale 
over a wide range of box sizes with definite exponents, 

m q (X) oc \ d + T (i\ (6) 

where A = fa/L is the quotient of box size fa to system 
size L. The exponents r(q) depend non-linearly on q and 
characterize the universality class of the system. Often 
the single exponent ao = dr(q) / dq(q = 0) is used in- 
stead of the entire r(q) spectrum, because it describes 
the scaling of the typical value of the squared amplitude 
with the system size, exp(ln l^l 2 ) cx L~ a ° . The results 
of numerical investigations of these exponent show that 
they do not depend on the size of W (i.e. W/E t , we will 
set E t — 1 in the following) U 1 However, we can observe 
that the range of box sizes for which scaling behavior is 
seen does change. For W < 1 the box probabilities scale 
over the entire range a < fa < L, where a is the lattice 
constant of the network. If we increase W, scaling dete- 
riorates for small box sizes, i.e. small A's. In Fig. 3 we 
show the ln[m 9 (Zf ) )]-ln A-plots for four different values of 
W at q = —0.5. We see that the deterioration of scaling 
shows a significant dependence on W indicating that the 
valid range of lengths for scaling follows a law of the kind 

a (J) < fa < L, (7) 

where a(x) is a function_we have to determine. 

If we take the valuestP of In A for each W at which 
the linear approximation becomes invalid and plot them 
against lnW, a linear dependence becomes apparent (see 
Fig. 4). A linear fit of the data yields a slope of X = 
1.36 ± 0.06. Consequently, the minimal length where 
scaling can be seen follows the following relation 



Another effect of the increase of W are the increasing 



sample to sample fluctuations and the rapid increase of 
the statistical error in the f(a) data, especially for large 
values of \q\. This all seems to be accompanied by a 
visible change of the characteristics of the wave function, 
as seen in Fig. 2. It is obvious that for W = 150 the wave 
function shows large areas, where the square amplitude 
changes only slightly, whereas for W = such areas of 
constant amplitude are much smaller. In fact one can 
see a steady increase of the size of these areas if one 
gradually increases W. 
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Fig. 4. The logarithms of the A's at which scaling deteriorates 
plotted against In W^. As a result for the linear fit we get a slope 
of 1.36 ±0.06 

If we take all these effects into account, we have to 
conclude that, apparently, the introduction of variable 
saddle point energies, i.e. W > 0, leads to a new length 
scale in our system. 

In order to understand how classical percolation can 
alter the shape of the quantum mechanical wave func- 
tions and influence their scaling behavior we will mpeat 
the discussion published previously in literature.Q'la'EJ' 
Take a look at those saddle points whose energy ut dif- 
fers more than E t from the electron energy E = 0. Of 
course, such saddle points occur only for W > E t . The 
transmission coefficients Tk at these saddle points are 
exponentially small and the amplitude of an incoming 
link is essentially transmitted to a single outgoing link 
(classical behavior). Therefore, on a path avoiding sad- 
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Fig. 5. Along the classical percolation path where the saddle 
point energies are large, > Et, the amplitude of the wave 
functions remain mainly constant. The resulting critical clusters 
form the links of the rcscaled network. 



die points with energy |ufc| < E t the amplitude of the 
wave function remains mainly constant. This leads to 
the formation of clusters where the wave function ampli- 
tude changes only slightly (see Fig. 5). The size of these 
clusters is obviously connected to classical percolation. 
The cluster correlation length £ p of clusters connected 
by saddle points whose energy is smaller than E t scales 
according to scaling theory (see e.g. refs. |l6|and [r7|) as 

This is the size of the largest clusters connected by those 
saddle points. If wc consider that for smaller clusters, 
which connect via these saddle points to other large or 
small clusters, the effect in terms of the change of wave 
function amplitude along the percolation path will be 
negligible, we see that the relevant cluster size is of order 
£ p . Accordingly, not all saddle points with Uf. < E t are 
relevant for the quantum mechanical effects, but only 
those critical saddle points that connect critical clusters 
and the scaling behavior of those clusters is given by 
(H). If we identify the clusters with the links of a new 
network and the saddle points with its nodes, we obtain a 
rescaled network model which has a renormalized lattice 
constant a 1 — a(W/E t ) and a new W — E t . In this way, 
we can treat network models with arbitrarily large finite 
W by jpapping it on the original Chalker-Coddington 
model Ja> Thus, we see that the influence of percolation 
is not on the overall scaling behavior but only on the 
range of the scaling. Its effects are mostly finite size 
effects which arise due to the decrease of the effective 
system size when we increase W . 

One very strong finite size effect is the shift of the 
critical energy E c when W is increased. This shift is not 
systematic but depends on the individual realization of 
the network. Due to the increasing range of the sad- 



dle point energies when W increases the average value 
of saddle points energy for a single system is no longer 
(u) = 0. Since the number of saddle points depends on 
the system size L, the number of random saddle point 
energies drawn from the square distribution may not be 
sufficient to guarantee that the deviation from the av- 
erage is small when the W approaches the system size 
causing the critical point to shift to E c = — (u) and con- 
sequently leading to strong fluctuations in the results for 
critical quantities. This suggests that any numerical in- 
vestigations should be confined to values of W that obey 
W < L. 

§4. Spectral Properties 

The next step in our investigation is to see what hap- 
pens to the spectral properties of the system when W is 
increased. In this case we are looking at the level spac- 
ing distribution P(s) and at the level number variance 
5^2 (N) = {(n — (n)) 2 ) of energy intervals containing on 
average N = (n) levels. .— . 

At the critical point P(s) has a unique shapdla^ which 
lies between the Poissonian behavior in the localized regime, 
P(s) — exp(— s), and that of the Gaussian unitary en- 
semble (GUE)-like in the metallic regime ( _P(s) = f|s 2 exp ( 
given by random matrix theory (RMT^tij For small s 
the level repulsion leads to a GUE-like behavior (P(s) oc 
s 2 ), whereas for larger s the tail of the distribution be- 
haves like P(s) oc exp (— ks), where k < 1. This mix- 
ture between GUE and Poisson is due to the multifractal 
structure of the wave functions which are neither homo- 
geneously smeared out over the entire system, as in the 
metallic regime, nor strongly localized. The direct con- 
nection between level distribution and_cpultifractal wave 
function was found by Chalker et al.t§> when they de- 
rived for the compressibility \ — lhrijv— >oc limL^oo dYi2 / dN 
at the critical point that \ — d ~2d 2 ' wncrc ^ i s the spatial 
and D 2 the fractal (correlation) dimension of the wave 
function. This equation could be vHrificd numerically for 
the Chalker-Coddington network^ and shows that we 
are again between the metallic (x = 0) and the localized 
(x = 1) behavior. If our arguments concerning the influ- 
ence of classical percolation are valid, we should obtain 
the same kind of results for the critical spectra regardless 
of the value of W. On the other hand, we have seen that 
increasing W causes finite size effects to appear faster 
which should also be visible in the spectral properties. 

We start by taking a look at the P(s) distribution for 
different values of W . In Fig. 6 we can see P(s) for W = 
to W — 40. In this case we chose the system size L = 80, 
but we obtain the same P(s) for other system sizes. The 
distribution is still system size independent at least for 
the range of system sizes we are able to investigate. This 
rules out finite size effects which could be expected due 
to the decrease of effective system size observed for the 
wave functions. Nevertheless, the distribution seems to 
become more and more Poissonian for larger W. This 
effect has already been seen by Ono et al.EP who used the 
Landau model for the quantum Hall effect. In this case 
an increase of the potential's correlation length, which is 
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Fig. 6. The P(s) spectrum for W = 0, 10, 20, 30, 40. The smooth 
curves show the Poisson and GUE distributions. The inset shows 
the same on a logarithmic scale. 



according to eq. (g) comparable to an increase of W, led 
to a similar result. In their publication Ono et al. fitted 
the P(s) distribution to the formula: 

As 13 exp(-Bs Q ), 



P(a) 



(10) 



where A and B are determined by the normalization con- 
ditions (s) = (1) = 1 and a and f} are the fitting param- 
eters. This way of—fitting was prompted by the results 
of Kravtsov et al.JiSEJ which predicted that /3 = 2 as 
for the Gaussian unitary ensemble (GUE), where P(s) — 
f|s 2 exp(-4/7rs 2 ), and a = 1+1/ (dv) « 1.21. Although, 
the numerical values for a do not fit that predictionEx 13 
we can still use eq. (|l(]) to compare our data with that 
in ref. ^l|. The result of such a fit to the distributions 
in Fig. 6 is shown in Table I. We can see that W w 30 
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P 







1.54 


±0.02 


2.03 


±0.03 


10 


1.71 


±0.04 


1.96 


±0.05 


20 


1.79 


±0.02 


1.30 


±0.02 


30 


1.94 


±0.03 


0.85 


±0.01 


40 


2.04 


±0.03 


0.59 


±0.01 


lc 


1.52 


±1.3 x 10" 4 


1.95 


±7.3 x 10" 4 


2l 


1.98 


±5.7 x 10~ 5 


1.03 


±1.6 x 10~ 3 



Table I. The results for the fit parameters a and (3 for different 
values of W. The last two entries give the values for different 
potential ranges determined by Ono et al. 



corresponds to the results obtained by Ono et al. for 
a potential correlation length of 21 c , whereas a correla- 
tion length of l c did not produce significant changes from 
the uncorrelated results (i.e. 14^ = 0). We would expect 
that, because the relevant microscopic length in the Lan- 
dau model is l c so that the potential correlation length 
has to be larger to show an appreciable effect. This cor- 
responds to the fact that W ~ E t also has no effect since 
the tunnel energy E t sets the energy scale for our system 
and W has to be larger in order to obtain classical saddle 



points. 

For both models it is very clear that the level repulsion 
for small values of s decreases with increasing W and 
correlation length, respectively. This is indicated by the 
increase of P(s = 0) and the decrease of the fit parameter 
P. 

In ref. [l3] it was found that the tail of P(s) shows an 
exponential decay of the form P-L s ) K cxp(— ks), which 
was predicted by Altshuler et alEJ* with k — 1/(2%) and 
confirmed for QHE systems by the numerical results in 
refs. [ll] and HJ. In Fig. 7 we show the behavior of the 
tail region of P(s) for different W. Although the tail 
seems to remain an exponential the factor k decreases 
with increasing W getting closer to the Poissonian be- 
havior, where k = 1. This would indicate that the com- 
1 
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Fig. 7. The tail of the P(s) spectrum for W = 0, 10, 20, 30, 40. 
The smooth curves are the linear fits for respective curves. 



pressibility \ also changes with W . In Fig. 8 we show 
the level number variance £2^) from which we obtain 
X by determining its slope. We can see that the slope 
increases with increasing W which would correspond to 
the decrease of k, but a look at Table II shows that the 
relation k — l/(2\) no longer holds in the case W > 0. 
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4.09 


±0.10 


0.127 


±0.002 


3.94 


±0.06 


10 


3.72 


±0.06 


0.190 


±0.003 


2.63 


±0.04 


20 


3.26 


±0.05 


0.254 


±0.001 


1.96 


±0.01 


30 


2.83 


±0.04 


0.291 


±0.002 


1.71 


±0.01 


40 


2.68 


±0.04 


0.313 


±0.002 


1.59 


±0.01 



Table II. The results for k and \ f° r different values of W . 



This is also true for the relation x 



d-D 2 
2d 



connect! 



the compressibility x to the multifractal exponent 
From the investigation of the critical wave functions we 
know that D2 does not change with W . 

Finally, observing that increasing W has a strong ef- 
fect on the shape of P(s) the question arises whether a 
function like 

1 



J (L,E) = -(s 2 ) 



L,E 1 



(11) 
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Fig. 8. The function E(JV) for different values of W. 



which depends on as well as describes the shape of P(s) 
still shows the same one-parameter scaling for W > as 
for W = 0. The insets in Fig. 9 shows Jo for W = 5, 15 
as a function of E for different system sizes. At the criti- 
cal point E c = all curves come together at a single size 
independent point which varies with W reflecting the de- 
pendence of the shape of P(s) on W. If one-parameter 
scaling still holds then a scaling function / of the follow- 
ing form should exist: 



J (E,L) = f(L/0, Z = Zo\E-E c \ 



(12) 



We determined f(x) by fitting a fourth degree polyno- 
mial in x — m(L/£o) + vlnE to our data using £o and 
v as fitting parameters. The values we found for v and 
the critical values of Jo are given in Table III. Although 
the critical values Jq (W) vary systematically with W the 
scaling exponent v shows no systematic dependence on 
W and is consistent. p4th, preidous results for v obtained 
by various studiesESEIIBESEP within the error bars. 
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2.1 


±0.3 


0.600 


±0.005 


2 


2.1 


±0.2 


0.601 


±0.01 


5 


2.2 


±0.1 


0.604 


±0.01 


10 


2.1 


±0.3 


0.611 


±0.01 


15 


2.3 


±0.1 


0.623 


±0.01 


30 


2.3 


±0.4 


0.664 


±0.01 


40 


2.2 


±0.3 


0.685 


±0.01 



Table III. The results for the scaling exponent v and the critical 
value of Jo for different values of W. 



§5. Conclusion 

In this paper we have discussed the influence of perco- 
lation effects on critical wave functions and the critical 
quasienergy spectrum of the Chalker-Coddington net- 
work model. Classical percolation effects appear in the 
generalized version of the model when the range W of 
the saddle point energies, i.e. node energies of the net- 
work, is increased from W = in the original model to 
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Fig. 9. The one-branch scaling functions Jo(-L/£) for four dif- 
ferent values of W resulting from the rescaling of the function 
Jq(E,L) shown in the insets. 



W > 1. 

We found that multifractal analysis of the critical wave 
functions results in the same exponents regardless of the 
value of W . Nevertheless, we can see the effect of clas- 
sical percolation in the appearance of a new microscopic 
length scale which increases with W and scales with an 
exponent X = 1.36 ± 0.06 consistent with the exponent 
Up = 4/3 of the classical correlation length. This con- 
firms the connection of the correlation length to that new 
length scale, which sets the minimum length for scaling 
and therefore reduces the effective size of the system, 
thereby causing increasing fluctuations when W grows. 

In the case of the spectral statistics of the quasienergies 
of the network the influence of an increasing W is more 
profound. An increasing shift from the critical statis- 
tics towards Poissonian behavior is observed for all the 
spectral quantities we have investigated. Although this 
shift does not seem to have an influence on the scaling 
exponent v it is clearly visible. Since the effect is not 
changing with system size, it can not be ruled as a mere 
finite size effect resulting from the decrease of effective 
system size observed for the critical wave functions. Fur- 
thermore, we observed that the changes that come with 
an increase of W are very similar to those connected with 
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the increase oiithe potential correlation length observed 
by Ono et al.B> 

It seems, although the multifractal behavior of the crit- 
ical wave functions is well understood, that more elab- 
orate analytical studies are indicated to understand the 
level statistics at the LD transition. 
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Appendix: Degeneracy of Eigenvalues 

The chiral structure of the Chalker-Coddington net- 
work leads to a pseudo degeneracy of the eigenvalues 
of its network operator U. For each eigenvector ty a with 
eigenvalue e m ° there exits an eigenvector ty a with eigen- 
value e ^= +7r ). 

In order to prove this, let us first point out the fact 
that the operator U maps vertical links only to hori- 
zontal links and vice versa. This means that U 2 maps 
vertical to vertical and horizontal to horizontal links, 
thereby creating two orthogonal £/ 2 -invariant subspaces. 
We can therefore define the projection operators P v and 
Ph which will project onto the two subspaces, respec- 
tively. 

Let \& a be an eigenvector of U with eigenvalue e lUa . 
Then ^f a is also an eigenvector of U 2 with eigenvalue 
e 2iu a -^Tg can wr j^ e \g a as a ii near combination of the 

two projections ^ v a = P v ^ a and = P h ^ a : 



ue„ 



he" 



(A-l) 



where v 



KJ/' 



„, = |**|, el = v-^l and e h a = 
h Since ^ v a and Vf^ are orthogonal to each other 

they are both also eigenvectors of U 2 with eigenvalue 
e 2iuj a _ ^his m eans that this eigenvalue is degenerate. 
Let us now construct the following eigenvector of U 2 : 



hel 



ve' 



^ Q is obviously orthogonal to $ tt . In order to see what 
will happen if we use U on ^q,, let us first use it on \& Q 
and keep in mind that Ue^ = e lulh e v a and Ue v a — e^^e 
where e 2 ^ = e 2lu]h = e 2lu <* . 

J* < he lUh ~ v - 



a » 



ve^e" 



(A-2) 



From eqs. flA2| ) and ( |Al| ) it follows that v = he i( - Ua ~ u ^ 

and h = ue'^°""' h ' and consequently 

UV a = he luJ "e'l~ve luJh e v a 

_ e 2i( Ua -u v ) he iu Ve h _ e 2i{ Ua -w h ) ve iu he v 



-e™° [-ve h c 



■he*] 



This means ^ a is an eigenvector of U with eigenvalue 



The same kind of argument leads to a further pseudo- 
degeneracy when we impose double periodic boundary 
conditions (torus) with system size L a multiple of 2. In 
that case each eigenphase uj a is accompanied by three 
other eigenphases w Q +7r/2, a; a +7r and uj a + 3n/2. This 
is caused by the chiral symmetry above and the bound- 
ary conditions which lead to the formation of two more 
(in this case U ) invariant subspaces within each of the 
J7 2 -invariant subspaces of vertical and horizontal lines. 
These subspaces consist of the sets combining every sec- 
ond vertical row or every second horizontal column, re- 
spectively Without the boundary conditions those sub- 
spaces mix. 
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